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Amendments to the Claims: 

This listing of claims will replace all prior versions, and listings, of claims in the application: 
Listing of Claims: 

Claim 1 (currently amended): A method for increasing the speed of line-by-line 
calculations of multiple overlapping absorption lines, said method comprising: 

(1) a step of dividing a domain of a Voigt function representing the shape of 
an absorption line into a first range around the peak of the Voigt function and a skirt 
portion not contained in the first range, replacing the first range with a cubic function 
whose function values and derivative values match with those of said Voigt function at 
the connection points with said Voigt function, and calculating the values and derivative 
values of said cubic function and the Voigt function in the skirt portion for each of first 
predetermined intervals; 

(2) a step of adding together the results of step (1 ) for a plurality of absorption 

lines; 

(3) a step of interpolating the resuUs of step (2) to calculate the function 
values and derivative values in second predetermined intervals smaller than said first 
predetermined intervals; 

(4) a step of dividing said first range into a second range near the peak and a 
skirt portion not contained in the second range, replacing said second range of a "function 
representing the difference between the Voigt function and said cubic function" with a 
second cubic function whose function values and derivative values match with those of 
said "function representing the difference between the Voigt function and said cubic 
function" at the connection points with said "function representing the difference between 
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the Voigt function and said cubic function", and calculating the values and derivative 
values of said second cubic function replacing said second range and the "function 
representing the difference between the Voigt function and sai d second cubic function" 
in the skirt portion for said second predetermined intervals; 

(5) a step of adding the results of step (4) to the results of step (3) for a 
plurality of absorption lines; 

(6) a step of interpolating the results of step (5) to calculate the function 
values and derivative values for intervals smaller than said second predetermined 
intervals; and 

(7) a step of adding the value s and derivative values of the "function 
representing the difference between the "" function representing the difference between 
the Voigt function and sai d second cubic function" " and said second cubic function" to 
the results of step (6) for a plurality of absorption lines in said second range. 

Claim 2 (original): A method in accordance with claim 1, wherein (8) the step 
of calculating function values and derivative values by interpolation is repeated while 
narrowing the interval until it becomes an interval of a minimum unit. 

Claim 3 (original): A method in accordance with claim 1 , wherein said steps (4) 
through (6) are repeated one or more times imtil the "intervals smaller than said second 
predetermined intervals" in step (6) reach third predetermined intervals. 

Claim 4 (original): A method in accordance with claim 1, wherein said 
predetermined interval is determined by using the following equation. 
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Said first predetermined intervals are f^'^dv. Here, j is a single-digit natural 
number, dv is the increment in wave number, and kmax is the largest natural number 
satisfying the relationship f^'^^^pdv < Vmax, where Vmax represents the maximum 
calculation range from the center of the absorption line, and p is 3, natural number for 
controlling the calculation precision. 



Claim 5 (original): A method in accordance with claim 3, wherein said 
predetermined interval is determined using the following equation. 

The third predetermined intervals are j^'"dv. Here, j is a single-digit natural 
number, dv is the increment of the wave number, and kmin is the maximum non-negative 
decimal number satisfying the relationship j^^pdv ^ a (a being approximately y/4) where 
Y is an approximate value of the full-width at half-maximum of the absorption line, and 
/? is a natural number for controlling the calculation precision. 



Claim 6 (original): A method in accordance with claim 1, wherein said second 
predetermined interval is determined by using the following equation. 

For the sub-function with the (^ - kmin + l)-th smallest width, the predetermined 
interval is/dv. Here,y is a single-digit natural number, dv is the increase in wave number 
and k is such that kmin ^ k< kmax. kmax is a natural number satisfying the relationship 
f^°^^^pdv = Fmox, where Vmax represents the maximum calculation range from the 
center of the absorption line, and is a natural number for controlling the calculation 
precision, kmin is the maximum non-negative decimal number satisfying the relationship 
j^'^pdv = a (a being approximately y/4), where y is an approximate value of the 
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full-width at half-maximum of the absorption line, and p is a natural nximber for 
controlling the calculation precision. 



Claim 7 (original): A method in accordance with any one of claims 1-3, wherein 
said interpolation is calculated by quartering the interpolation interval, with the function 
values represented hyy^, y^^ and the derivative values represented by 3/^', y^\ y^* at the 
points where the interpolation interval (jc^ Xy) is divided, using the function values j/^ y^ 
and function derivative values y^, y/ at Xq, X;, with the below-given Equation (1) as a 
function value interpolation equation, Equation (2) as a function derivative value 
interpolation equation and c as a non-negative decimal fraction. 



yb 

\yc) 



j_ 

64 



^54-6£r 10 + 6f 9(1-.f) -3(1-^^ 

32 32 8(1 -f) -8(1- 
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-24 + 8f 24-8^ -4(l-f) -4(l-f) 



^-18H-2f 18-26: -5(1-^) 3(l-£-) 



yo 
yi 



(2) 



Claim 8 (original): A method in accordance with any one of claims 1 -3, wherein 
said interpolation is calculated by dividing the interpolation interval into five parts, with 
the function values represented by^„, y^^ y^ and the derivative values represented by 
ya' yb' yc- yd at the points where the interpolation interval {xg, x,) is divided, using the 
function values y, and function derivative values 3^0'. yi'^^Xg, x„ with the below-given 
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Equation (3) as a function value interpolation equation. Equation (4) as a function 
derivative value interpolation equation and e as a non-negative decimal fraction. 
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(4) 



Claim 9 (original): A method in accordance with any one of claims 1-8, for 
increasing the speed of line-by-line calculations, when assuming the Voigt function to be 
K(x, y) and the shape of the absorption line displaced from the Voigt profile of the 
absorption line to be K(x. y) + f(x), by replacing 

i^ix^y) 

with: 

K{x,y)= AK{x,yyBf{x) 

and 

dK{x,y) 
dx 
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with 

dK(x,y) _^ dK{x,y) ^^dfjx) 
dx dx dx 

Claim 10 (original): A method in accordance with any one of claims 1-8, for 
increasing the speed of line-by-line calculations, when assuming the Voigt function to be 
K(x, y) and the shape of the absorption line displaced from the Voigt profile to be K(x, 
y)f(x)^ by replacing 

K{x.y) 

with: 

K{x,y)= K{x,y)f{x) 

and 

dK{x,y) 
dx 

with 

dx dx dx 

Claim 11 (original): A method in claim 10, for increasing the speed of 
line-by-line calculations, by using, in a sub-Lorentzian correction, 

K{x, y)=K (x, y)A exp (- B\x\) 
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and 



= ^ exp(- K{x, yl- sgn {x)AB exp(- 

Claim 12 (original): A method in accordance with claim 9, for increasing the 
speed of line-by-line calculations, by replacing, in a line-mixing correction. 



with: 



K{x,y)=AK{x,y)+BL{x,y) 



and 



dK{x,y) 
dx 



with 



dKix,y) 
dx 



= -2 



(Ax + By)K{x, y)-{Ay - Bx)l{x, y)- 



B 



Here, L(x, y) is the imaginary component of the function w(z) (the real part is the Voigt 
function), where complex number z = x + iy, defined by the following equation. 

{erfc(z) is a complex complementary error function). 
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Claim 1 3 (currently amended): A program for increasing the speed of line-by-line 
calculations of multiple overlapping absorption lines, said program comprising: 

(1) a step of dividing a domain of a Voigt function representing the shape of 
an absorption line into a first range around the peak of the Voigt function and a skirt 
portion not contained in the first range, replacing the first range with a cubic function 
whose function values and derivative values match with those of said Voigt function at 
the connection points with said Voigt function, and calculating the values and derivative 
values of said cubic function and the Voigt function in the skirt portion for each of first 
predetermined intervals; 

(2) a step of adding together the results of step ( 1 ) for a plxirality of absorption 

lines; 

(3) a step of interpolating the results of step (2) to calculate the function 
values and derivative values in second predetermined intervals smaller than said first 
predetermined intervals; 

(4) a step of dividing said first range into a second range near the peak and a 
skirt portion not contained in the second range, replacing said second range of a "function 
representing the difference between the Voigt function and said cubic function" with a 
second cubic function whose function values and derivative values match with those of 
said "function representing the difference between the Voigt function and said cubic 
function" at the coimection points with said "function representing the difference between 
the Voigt function and said cubic function", and calculating the values and derivative 
values of said second cubic function replacing said second range and the "function 
representing the difference between the Voigt function and sai d s e cond cubic function" 
in the skirt portion for said second predetermined intervals; 
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(5) a step of adding the results of step (4) to the results of step (3) for a 
plurality of absorption lines; 

(6) a step of interpolating the results of step (5) to calculate the function 
values and derivative values for intervals smaller than said second predetermined 
intervals; and 

(7) a step of adding the value s and derivative values of the "function 
representing the difference between the " "function representing the difference between 
the Voigt function and sai d sec o nd cubic function" " and said second cubic function" to 
the results of step (6) for a plurality of absorption lines in said second range. 

Claim 14 (original): A program in accordance with claim 13, wherein the step 
of calculating function values and derivative values by interpolation is repeated while 
narrowing the interval imtil it becomes an interval of a minimum unit. 

Claim 15 (original): A program in accordance with claim 13, wherein said steps 
(4) through (6) are repeated one or more times until the "intervals smaller than said 
second predetermined intervals" in step (6) reach third predetermined intervals. 
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